Rubber friction and tire dynamics 
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We propose a simple rubber friction law, which can be used, e.g., in models of tire (and vehicle) 
dynamics. The friction law is tested by comparing numerical results to the full rubber friction theory 
(B.N.J. Persson, J. Phys.: Condensed Matter 18, 7789 (2006)). Good agreement is found between 
the two theories. 

We describe a two-dimensional (2D) tire model which combines the rubber friction model with 
a simple mass-spring description of the tire body. The tire model is very flexible and can be 
used to calculate accurate /x-slip (and the self-aligning torque) curves for braking and cornering or 
combined motion (e.g., braking during cornering). We present numerical results which illustrate the 
theory. Simulations of Anti-Blocking System (ABS) braking are performed using two simple control 
algorithms. 



1. Introduction 

Rubber friction is a topic of huge practical importance, 
e.g ., for tires, rubber seals, conveyor belts and syringes[l|- 
17[. In most theoretical studies rubber friction is de- 
scribed using very simple phenomenological models, e.g., 
the Coulombs friction lavif with a friction coefficient which 
may depend on the local sliding velocity. However, as we 
have shown earlier 6j, rubber friction depends on the his- 
tory of the sliding motion (memory effects), which we 
have found to be crucial for an accurate description of 
rubber friction. For rubber sliding on a hard rough sub- 
strate, the history dependence of the friction is mainly 
due to frictional heating in the rubber-substrate contact 
regions. Many experimental observations, such as an ap- 
parent dependence of the rubber friction on the normal 
stress, can be attributed to the influence of frictional 
heating on the rubber friction. 

A huge number of papers have been published related 
to tire dynamics, in particular in the context of Anti- 
Blocking System (ABS) braking models. The "heart" in 
tire dynamics is the road-rubber tire friction. Thus, un- 
less this friction is accurately described, no tire model, 
independent of how detailed the description of the tire 
body may be, will provide an accurate picture of tire dy- 
namics. However, most treatments account for the road- 
tire friction in a very approximate way. Thus, many "ad- 
vanced" finite element studies for tire dynamics account 
for the friction only via a static and a kinetic rubber fric- 
tion coefficients. In other studies the dynamics of the the 
whole tire is described using interpolation formulas, e.g., 
the "Magic Formula" {4J , but this approach require a very 
large set of measured tire properties (which are expensive 
and time-consuming to obtain), and cannot describe the 
infiuence of history (or memory) effects on tire dynamics. 

In this paper we first propose a very simple rubber fric- 
tion law (with memory effects) which gives nearly iden- 
tical results to the full model developed in Ref. |6[. We 
also develop a 2D-tire model which combines the rubber- 
road friction theory (which accounts for the flash temper- 
ature) with a simple two-dimensional (2D) description of 



the tire body. We believe that the most important aspect 
of the tire body is its distributed mass and elasticity, and 
this is fully accounted for in our model. One advantage of 
the 2D-model over a full 3D-model is that one can easily 
impose any foot-print pressure distribution one like (e.g., 
measured pressure distributions), while in a 3D-model 
the pressure distribution is fixed by the model itself. This 
allow a detailed study on how sensitive the tire dynamics 
depends on the nature of the foot print pressure distribu- 
tion. The tire model is illustrated by calculating /i-slip 
curves, and with simulations of ABS braking using two 
different control algorithms. 

2. Rubber friction 

Rubber friction depends on the history of the sliding 
motion. This is mainly due to the flash temperature: the 
temperature in the rubber-road asperity contact regions 
at time t depends on the sliding history for all earlier 
times t' < t. This memory effect is crucial for an ac- 
curate description of rubber friction. We illustrate this 
effect in Fig. [T] for a rubber tread block sliding on an 
asphalt road surface. We show the (calculated) kinetic 
friction coefficient for stationary sliding without (blue 
curve) and including the fiash temperature (red curve) 
as a function of the velocity v of the bottom surface of 
the rubber block. The black curves shows the effective 
friction during non-stationary sliding experienced by a 
rubber tread block during braking at various slips (slip 
values from 0.005 to 0.09). Note that because some finite 
sliding distance is necessary in order to fully develop the 
flash temperature, the friction acting on the tread block 
initially follows the blue curve corresponding to "cold- 
rubber" (i.e., negligible flash temperature). Thus, it is 
not possible to accurately describe rubber friction with 
just a static and a kinetic friction coefficients (as is often 
done even in advanced tire dynamics computer simula- 
tion codes) or even with a function fi(v) which depends 
on the instantaneous sliding velocity v{t). Instead, the 
friction depends on u(t') for all times t' < t. 

As a background to what follows, we first review the 
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FIG. 1: Red and blue lines: the kinetic friction coefBcient 
(stationary sliding) as a function of the logarithm (with 10 
as basis) of the sliding velocity. The blue line denoted "cold" 
is without the flash temperature while the red line denoted 
"hot" is with the flash temperature. Black curves: the effec- 
tive friction experienced by a tread block as it goes through 
the foot print. For the car velocity 27 m/s and for several 
slip values 0.005, 0.0075, 0.01, 0.03, 0.05, 0.07 and 0.09. Note 
that the friction experienced by the tread block first follows 
the "cold" rubber branch of the steady state kinetic friction 
coefficient and then, when the block has slip a distance of 
order the diameter of the macroasperity contact region, it 
follows the "hot" rubber branch. 



rubber friction theory (see [3,Q for details). It is assumed 
that all energy dissipation arises from the viscoelastic de- 
formations of the rubber surface by the road asperities. 
An asperity contact region with the diameter d gives rise 
to time dependent (pulsating) deformations of the rub- 
ber which is characterized by the frequency lo — v/d, 
where v is the sliding velocity. The viscoelastic deforma- 
tion (and most of the energy dissipation) extend into the 
rubber by the typical distance d. So most of the energy 
dissipation occur in a volume element of order d^ . In or- 
der to have a large asperity-induced contribution to the 
friction the frequency uj should be close to the maximum 
of the tanS = linE{Lu) /KeE{Lu) curve. Here E{uj) is the 
viscoelastic modulus of the rubber. In reality there will 
be a wide distribution of asperity contact sizes, so there 
will be a wide range of perturbing frequencies, say from 
Wo to uji, see Fig. [51 A large friction requires that tan^ is 
as large as possible for all these perturbing frequencies. 

The temperature dependence of the viscoelastic mod- 
ulus of rubber-like materials is usually very strong, and 
an increase in the temperature by 10 °C may shift the 
tan(5 curve to higher frequencies by one frequency decade. 
This will usually reduce the rubber friction, see Fig. [51 

Real surfaces have a wide distribution of asperity sizes. 
The best picture of a rough surface is to think about it 
as big asperities on top of which occur smaller asperi- 
ties on top of which occur even smaller asperities ... . 
This is illustrated in Fig. [31 for a case where roughness 
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FIG. 2: When the temperature increases the tanS — 
ImE /KeE spectra shift to higher frequencies, which result 
in a decrease in the rubber friction. We assume the road as- 
perities gives rise to pulsating frequencies in the range ujq and 

ci;i. 



rubber 



5.maj!..ertergy... ..y 



*■ .|!„.fyi-^sipktiih 

^li 11 m"'M., 11 



lirgi'eriirgf V 

i|/iisiiati:h m \ 




FIG. 3: The dissipated energy per unit volume is highest in 
the smallest asperity contact regions. 



occur on two length scales. To get the total energy dis- 
sipation during sliding on a real surface one need to sum 
up the contribution from asperity induced deformation 
of the rubber on all (relevant) length scales. It is impor- 
tant to note that all length scales are a priory equally 
important [5|. 

Temperature has a crucial influence on rubber friction. 
The viscoelastic energy dissipation, which is the origin of 
the rubber friction in my model, result in local heating 
of the rubber in exactly the region where the energy dis- 
sipation occur. This result in a temperature increase, 
which becomes larger as we observe smaller and smaller 
asperity contact regions. This local (in time and space) 
temperature increase, resulting from the local viscoelas- 
tic energy dissipation, is referred to as the flash temper- 
ature. The flash temperature has an extremely impor- 
tant influence on the rubber friction. This is illustrate 
in Fig. [1] where we show the (calculated) steady state 
kinetic friction coefhcient when a block of tread rubber 
is sliding on an asphalt road surface. The upper curve is 
the result without accounting for the flash temperature, 
i.e., the temperature is assumed to be equal to the back- 
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FIG. 4: The contact region between a tire an a road sur- 
face. At low magnification (" < 1 it appears that the tire is in 
complete contact with the road but as the magnification in- 
creases, the contact area continuously decreases as indicated 
in the figure. 



ground temperature Tq everywhere. The lower curve is 
the result including the flash temperature. Note that for 
sliding velocities v > 0.001 m/s the flash temperature 
result in a lowering of the sliding friction. For velocities 
V < 0.001 m/s the produced heat has enough time to 
diffuse away from the asperity contact regions, and the 
flash temperature effect is negligible. 

In the rubber friction theory the concept of the 
macroasperity contact region is of crucial importance. 
Let us study the footprint contact region between a tire 
and a road surface at different magnification C. At low 
magnification the road surface appears smooth and the 
contact between the tire and the road appears to be com- 
plete within the footprint area as in Fig. |4l[a). However, 
when we increase the magnification ^ we start to observe 
non-contact regions as in Fig. |H[b) . At high enough mag- 
nification we observe isolated contact regions as in Fig. 
|4l^c), which, when the magnification increases even fur- 
ther, break up into even smaller contact regions as in Fig. 
m^d). We denote the contact regions observed in Fig. Hfc) 
as the macroasperity contact regions (with the avera ge d i- 
ameter D) (the exact definition is given in Ref. |6l.ll8|). 
If the nominal pressure in the tire-road contact region 
is small the macroasperity contact regions will be well 
separated, but the separation between the microasperity 
contact regions within the macroasperity contact regions 
is in general very small. When calculating the flash tem- 
perature effect we have therefore smeared out the heat 
produced by the microasperity contact regions uniformly 
within the macroasperity contact regions. Typically for 
road surfaces D « 0.1 — 1 cm, and the fraction of the 
tread block surface occupied by the macroasperity con- 
tact regions are typically between 10% and 30%. 



In the friction theory developed in Ref. [6| only the 
surface roughness with wave vectors q < qi are assumed 
to contribute to the friction. For clean road surfaces 
we determine the cut-off wavevector qi by a yield con- 
dition: We assume that the local stress and temperature 
in the asperity contact regions on the length scale 1/gi 
are so high that the rubber bonds break resulting in a 
thin modified (dead) layer of rubber at the surface re- 
gion of thickness w 1/gi. From this follows the following 
observations: 

(a) The rubber friction on clean road surfaces after 
run-in is rather insensitive to the road surface. This has 
been observed in several series of experimental studies 
(not shown), and is in accordance with the present the- 
ory. This can be understood as follows. The cut-off qi on 
surfaces with smoother, less sharp roughness or surfaces 
where the roughness occur at shorter length scales will be 
larger (i.e. the cut-off wavelength Ai = 27r/gi smaller) 
than for road surfaces with larger roughness in such a 
way as the stress and temperature increase in the asperity 
contact regions which can be observed at the resolution 
Ai (or magnification ( = qi/qo) are roughly the same on 
all surfaces. This implies that a larger range of rough- 
ness will contribute to the rubber friction on "smoother" 
surfaces as compared to more rough surfaces. As a con- 
sequence, the friction (after run-in) may vary very little 
between different (clean) road surfaces. 

(b) On contaminated road surfaces, the cut-off qi may 
be determined by the nature of the contamination. In 
this case, if the cut-off is fixed (e.g., determined by, say, 
the size of contamination particles) one may expect much 
larger variation in the friction coefficient between differ- 
ent road surfaces, and also a larger variation between 
tired with different types of tread rubber. 

3. Phenomenological rubber friction law 

In tire applications, for slip of order 5 — 10% and typical 
footprint length of order 10 cm, the slip distance of a 
tread rubber block in the footprint will be of order 1 cm, 
which typically is of order the diameters D of the macro 
asperity contact regions. As discussed above, as long as 
the slip distance s{t) is small compared to D one follow 
the cold rubber branch of the steady state relation ij,{v) 
so that /i(t) w A'coid(i'(i)) for the slip distance s{t) << 
D. When the tread block moves towards the end of the 
footprint the slip distance s{t) may be of order of (or 
larger than) D, and the friction will follow the hot branch 
of the /i(w) relation i.e., /i(t) « /ihot(w(0) fo^' •s(^) > D. 
We have found that the following (history dependent) 
friction law gives nearly the same result as the full theory 
presented in Ref. [6|: 
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FIG. 5: The frictional shear stress acting on a tread block as 
a function of time for many slip values: 0.005, 0.0075, 0.01, 
0.03, 0.05, 0.07, 0.09, 0.12, 0.15 and 0.25. For the car velocity 
27 m/s and tire background temperature To = 60 °C. For the 
ID tire model using the full friction model (green curves) and 
the cold-hot friction law (1) (red curves). For a passenger car 
tread compound. 



where v(t) is the instantaneous sliding velocity and s(t) 
the sliding distance, and sq « 0.2D. We will refer to (1) 
as the cold-hot friction law. The length D depends on the 
rubber compound and the road surface but is typically in 
the range D f* 0.1 — 1 cm. Using the full friction theory 
one can easily calculate the functions ^J-coldiv{t),TQ) and 
IJ'hot{v{t),To) and the length D. 

To demonstrate the accuracy of the cold-hot rubber 
friction law (1), let us study the dynamics of one tread 
block as it pass through the tire-road footprint. In Fig. [5] 
we show the the frictional shear stress acting on a tread 
block as a function of time for many slip values: 0.005, 
0.0075, 0.01, 0.03, 0.05, 0.07, 0.09, 0.12, 0.15 and 0.25. 
Note that the cold- hot friction law (1) (red curves) gives 
nearly the same result as for the full friction model (green 
curves). In Fig. [S]we show the /x-slip curve. Again the 
cold-hot friction law (1) (red curve) gives nearly the same 
result as the full friction model (green curve). 

4. Tire dynamics 

All the calculations presented in this section have been 
obtained using the 2D tire model described below, with 
the tire body optimized using experimental data for a 
passenger car tire. The viscoelastic springs associated 
with this tire body are kept fixed in all the calculations. 
Thus the model calculations does not take into account 
the changes in the tire body viscoelastic properties due 
to variations in the tire (background) temperature, or 
variations in the tire inflation pressure (which affect the 
tension in the tire walls). In principle both effects can be 
relatively simply accounted for in the model, but have not 
been included so far. All the calculations presented below 
are obtained using the measured viscoelastic modulus of 
a passenger car tire tread compound, and for the Estoril 




FIG. 6: The /.i-slip curve for the ID tire model using the full 
friction model (green curve) and the cold-hot friction law (1) 
(red curve). For a passenger car tread compound. 
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FIG. 7: 2D-model of a tire. The car velocity Vc points in 
another direction than the rolling direction, giving a non-zero 
cornering angle 0. 



racer track. 

We use a 2D-description of the tire body as indicated 
in Fig. [7l We introduce a coordinate system with the 
y-axis in the transverse direction and the x-axis along 
the longitudinal (rolling) direction. We consider the tire- 
road system in a reference frame where the road is sta- 
tionary. The car velocity Vc, the rolling velocity vb. and 
the cornering angle determines the transverse Vy and 



longitudinal v^ slip velocities: 



Vx 



Vy — i;cSin6' 



UcC0s6' — Vb. 



The longitudinal slip is defined by 

Vx ticC0s6' — Ur 

s = — 

WcCOS^ WcCOs6' 

When the cornering angle 9 = this equation reduces to 
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Note that Vx and Vy are also the velocities of the tire rim, 
and that the footprint moves (in the rolling direction) 
with the velocity UcCos0 relative to the road, and with 
the rolling velocity WcCOS^ — Vx = wr relative to the rim. 

We describe the tire body as a set of mass points 
connected with viscoelastic springs (elasticity k and vis- 
cous damping 7). The springs have both elongation and 
bending elasticity, denoted by k and fc', respectively, and 
the corresponding viscous damping coefficients 7 and 7'. 
We assume Nx and Ny tire body blocks along the x- 
and y-directions and let x^j = [xij,yij) denote the dis- 
placement vector of tire body block (i,j) {i = 1, ■■■,Nx, 
j = l,...,Ny). Since the tire is a torus shaped object 
we must assume periodic boundary conditions in the x- 
direction so that xn^+i,j — xij and yM^+ij = yi,j- 

For stationary tire motion we have the following 
boundary conditions. For i = 0, ..., Nx + 1: 

ViO = Vyt, Xio = Vxt 



yi,Ny+l = Vyt, Xi^Ny+1 = Vxt 

FoTJ = l,...,Ny: 

VN^ + lJ = 2/lj> ^N^ + l,j = Xij 

If the mass of tire body element (i,j) is denoted by 
mj, we get for i = 1, ...,Nx, j = l,...,Ny: 

m-jjj.j = Fy,j + kyjiy^j^i^ - y^-) + kyj+i{y,^.i+i - y.-j) 



+Kjiyr+h3 +y»-ij -22/jj) 



+i'xjiyt+h3 +yt-ij -'^mj) 



FIG. 8: Three tire-road footprints used in the computer sim- 
ulations. The footprints are 0.2 m wide and the normal pres- 
sure in the footprints is constant at p = 0.1 MPa or 0.3 MPa. 
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In the equations above, Fxij and Fyij are the force com- 
ponents (in the x- and y- direct ions, respectively) acting 
on the tire body block (z, j) from the tread block (i,j). 
Thus, Fij = {Fxij, Fyij) is nonzero only when {i,j) is 
in the tire tread area. The tire body viscoelastic spring 
parameters (A:, 7) and (A;', 7') in the tread area and in 
the side wall area (8 parameters) have been optimized in 
order to reproduce a number of measured tire properties 
(e.g., the longitudinal and transverse tire stiffness val- 
ues for three tire loads, and the frequency and damping 
of the lowest longitudinal and transverse tire vibrational 
modes). The optimization has been performed usiiig the 
amoeba method of multidimensional minimization|l9|. 

4.1 Dependence of the /i-slip curve on the shape 
of the tire-road footprint 

Here we will study how the tire dynamics depends on 
the shape of the tire-road footprint. We consider the 
rectangular, elliptic and trapezoid footprints shown Fig. 
[5] We assume first that the pressure in the footprint is 
uniform and equal to p = 0.1 MPa, and that the tire load 
is equal to 2000 N in all cases. Thus all the foot prints 
have the same area. 

In Fig. [9] we show the /i-slip curves for the three foot- 
prints shown in Fig. [51 and in Fig. [TU] we show the 
corresponding /i-slip angle curves. It is remarkable how 




FIG. 9; The /i-slip curves for the three footprints shown in 
Fig. |8] For the rubber background temperature To = 80 °C 
and the car velocity 27 m/s. For Fn = 2000 N and p = 
0.1 MPa. 




FIG. 11: The self-aligning moment for the trapezoid footprint, 
as a function of the longitudinal slip. For the rubber back- 
ground temperature To = 80 °C and the car velocity 27 m/s. 
For Fn = 2000 N and p = 0.1 MPa. 
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FIG. 10: The /x-slip angle curves for the three footprints 
shown in Fig. |8l For the rubber background temperature 
To = 80 °C and the car velocity 27 m/s. For Fn = 2000 N 
andp = 0.1 MPa. 
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FIG. 12: The self-aligning moment for the three footprints 
shown in Fig. [8] as a function of the coornering angle. For 
the rubber background temperature To = 80 °C and the car 
velocity 27 m/s. For Fn = 2000 N and p = 0.1 MPa. 



insensitive the results are to the shape of the footprint. 
Thus we may state that ^-slip curves, and hence tire dy- 
namics, depends very weakly on the shape of the tire-road 
footprint, assuming everything else the same. 

In Fig. [TT] we show the the self-aligning moment for 
the trapezoid footprint profile in Fig. |S1 as a function 
of the longitudinal slip (for zero cornering angle) . Note 
that for the rectangular and elliptic footprint the self 
ahgning moment vanish (not shown). This is expected 
because of the mirror symmetry of the footprint in the 
X-axis through the center of the footprint. However, the 
trapezoid footprint does not exhibit this symmetry (see 
Fig. [5]), and the self aligning moment is non- vanishing in 
this case. 

In Fig. [T2] we show the self-aligning moment for the 
three footprints shown in Fig. [51 as a function of the 
cornering angle (for zero longitudinal slip). In this case 



the self aligning moment is non-vanishing in all cases. It 
is remarkable, however, how insensitive the results are to 
the shape of the footprint. 

For small 9 the self aligning moment is positive. This is 
due to the gradual (nearly linear) build up of the trans- 
verse stress from the inlet of the footprint to the exit. 
Thus the center of mass of the frictional stress distribu- 
tion is located closer to the exit of the footprint giving a 
positive self aligning moment. For large 6 the self align- 
ing moment is negative. In the present model this is due 
to the flash temperature effect: when a tread block enter 
the footprint the rubber is "cold" , and initially the rub- 
ber friction is high. After a short slip distance the full 
flash temperature is built up and the rubber friction is 
smaller. This will result in a frictional stress distribution 
which is larger close to the inlet. This in turn result in a 
negative self aligning moment. 

In the study above we have assumed a constant pres- 
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FIG. 13: The /i-slip angle curves for the three footprints 
shown in Fig. [8] and for the fongitudinal slip s = and s = 
0.03. For the rubber background temperature To = 80 °C and 
the car velocity 27 m/s. For _Fn = 2000 N and p = 0.1 MPa. 
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FIG. 14: The self-ahgning moment for the three footprints 
shown in Fig. [8] and for the longitudinal slip s = 0.03. For 
the rubber background temperature To = 80 °C and the car 
velocity 27 m/s. For _Fn = 2000 N and p = 0.1 MPa. 



sure in the footprints. In reality, the pressure will be 
slightly larger at the inlet than at the exit of contact 
with the road (this is the case also during pure rolling 
and is related to the rolling resistance). This asymme- 
try will give an additional (negative) contribution to self 
aligning moment for large slip. 

In Fig. [13] we show the /i-slip angle curves for the 
three footprints shown in Fig. |51 and for the longitudi- 
nal slip s — and s = 0.03. Note, in accordance with 
experimental observations, for the combined slip ^iy{0) 
is smaller than for the case when the longitudinal slip 
vanish. Again there is very little influence on the shape 
of the footprint. However, the self aligning moment will 
now depend strongly on the shape of the footprint. This 
is shown in Fig. [14] for the case s — 0.03. Note that 
for the trapezoid footprint the self aligning moment is 
much larger than for the elliptic and rectangular foot- 
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FIG. 15: The ^-slip angle curves (a), and the self aligning mo- 
ment (b) for the elliptic footprint for the tire load Fn = 3000, 
5000 and 7000 N, and the footprint pressure p = 0.3 MPa. 
For the rubber background temperature To = 80 °C and the 
car velocity 27 m/s. 



prints. This is due to the contribution from the longi- 
tudinal stress component which gives rise to a net longi- 
tudinal force centered to the right of the mid-line of the 
tire. 

Fig. [TSl shows the /i-slip angle curves (a), and the self 
aligning moment (b) for the elliptic footprint for the tire 
load F^ = 3000, 5000 and 7000 N, and the footprint pres- 
sure p = 0.3 MPa. Note that as the load increases the 
foot print becomes longer which result in a decrease in 
the maximum friction coefficient, which agree with exper- 
imental observations. This load-dependence is not due to 
an intrinsic pressure dependence of the rubber friction co- 
efficient (which was kept constant in our calculation) , but 
a kinetic effect related to the build up of the flash tem- 
perature in rubber road asperity contact regions during 
slip. To understand this in more detail, consider again 
Fig. [2 

The red and blue lines in Fig. [T] show the kinetic fric- 
tion coefficient (stationary sliding) as a function of the 
logarithm of the sliding velocity. The upper line denoted 
"cold" is without the flash temperature while the lower 
line denoted "hot" is with the flash temperature. The 
black curves show the effective friction experienced by 



a tread block as it goes through the footprint. Results 
are shown for several slip values 0.005, 0.0075, 0.01, 0.03, 
0.05, 0.07 and 0.09. Note that the friction experienced 
by the tread block first follows the "cold" rubber branch, 
and then, when the block has slipped a distance of order 
the diameter D of the macroasperity contact region, it 
follows the "hot" rubber branch. Based on this figure it is 
easy to understand why the maximum friction coefficient 
increases when the length of the footprint decreases: If 
Wsiip is the (average) slip velocity of the tread block, then 
in order to fully build up the fiash temperature the fol- 
lowing condition must be satisfied: WsiipisUp ~ D, where 
D is the diameter of the macroasperity contact region. 
Since the time the rubber block stay in the footprint 
^siip — L/vYi (where L is the length of the footprint and 
vb. the rolling velocity) we get WsUp « vb.{D/L). Thus, 
when the length L of the footprint decreases, the (aver- 
age) slip velocity of the tread block in the footprint can 
increase without the slip distance exceeding the diameter 
D of the macro asperity contact region. Thus at the slip 
corresponding to the maximum of the /i-slip curve, for 
a short footprint, as compared to a longer footprint, the 
tread block will follow the "cold" rubber branch of the 
(steady state) /^-slip curve to higher slip velocities before 
the fiash temperature is fully developed, resulting in a 
higher (maximal) tire-road friction for a short footprint 
as compared to a longer footprint. 

In Fig. [16] we show snapshot pictures of the tire body 
deformations for three cases, namely with the slip s and 
the coornering angle 9 given by {s,9) — (0.05,0) (left), 
(0.03,5°) (middle) and (0,5°) (right). The open squares 
denote the position of rubber elements of the undeformed 
tire body, and the filled squares underneath denote the 
position of the same tire body elements of the deformed 
tire. 

In Fig. [17] we show snapshot pictures of the tire body 
deformations for the rectangular, trapezoid and elliptic 
contact area with contact pressure 0.3 MPa and normal 
load Fn = 7000 N. In aU cases the slip s = 0.05 and 
the coornering angle 9 — 0. Note how insensitive the 
deformation field is to the shape of the footprint. This is 
caused by the high stiffness of the tire body in the tread 
area. 

4.2 Dependence of the /i-slip curve on the size 
of the tire-road footprint 

In Fig. [18] we show /i-slip curves for the rectangular 
footprint for the contact pressures p = 0.1, 0.125, 0.15 
and 0.3 MPa. The tire load is fixed at Fn = 2000 so the 
different contact pressures corresponding to the footprint 
length L = 10.2, 8.1, 6.8 and 3.4 cm, respectively. Note 
that increasing the tire footprint pressure decreases the 
length of the footprint, which decreases the tire longitu- 
dinal stiffness Cx (determined by the initial slope of the 
Ma;(s)-curves), and also the maximum of the /i-slip curves. 
In practice the pressure in the tire-road footprint can be 
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FIG. 16: Snapshot pictures of the tire body deformations for 
shp s and the coornering angle given by (s^O) = (0.05,0) 
(left), (0.03,5°) (middle) and (0,5°) (right). The open 
squares denote the position of rubber elements of the unde- 
formed tire body, and the filled squares underneath denote the 
position of the same tire body elements of the deformed tire. 
For the trapezoid contact area with contact pressure 0.1 MPa 
and normal load -Fn = 2000 N. For the rubber background 
temperature To = 80 °C and the car velocity 27 m/s. 



changed by changing the tire inflation pressure, but in 
this case one expect also a change in the tire body stiff- 
ness. In the tire model we use this effect is not included 
at present. Thus, at least for passenger car tires, when 
the inflation pressure increases the longitudinal tire stiff- 
ness Cx usually first increases and then, at high enough 
inflation pressure, decreases. This is consistent with Fig. 
1181 which shows initially a very small change in the tire 
stiffness as the footprint pressure p increases, so that for 
small (but not too small) p, the stiffening of the tire body 
may dominate over the contribution from the change in 
the footprint, so that initially Cx increases with increas- 
ing p. 

In Fig. [TO] we show snap shot pictures of the tire body 
deformations for three cases, namely for the external load 
Fn = 3000, 5000 and 7000 N. In all cases the slip s = 0.05 
and the coornering angle ^ = 0. The results are for the 
rectangular footprint with contact pressure 0.3 MPa. 

In Fig. [20] we show similar results as in Fig. \W\ but 
now for the elliptic contact area with contact pressure 
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FIG. 17: Snapshot pictures of the tire body deformations for 
the rectangular, trapezoid and elhptic contact area with con- 
tact pressure 0.3 MPa and normal load Fn = 7000 N. For 
the slip s = 0.05 and the coornering angle ^ = 0. The open 
squares denote the position of rubber elements of the unde- 
formed tire body, and the filled squares underneath denote 
the position of the same tire body elements of the deformed 
tire. For the rubber background temperature To = 80 °C and 
the car velocity 27 m/s. 




FIG. 18: The /i-slip curves for the rectangular footprint shown 
in Fig. [S] For the contact pressures p = 0.1, 0.125, 0.15 
and 0.3 MPa corresponding to the footprint length L = 10.2, 
8.1, 6.8 and 3.4 cm. For the rubber background temperature 
To = 80 °C, the tire load Fn = 2000, and the car velocity 
27 m/s. 



FIG. 19: Snapshot pictures of the tire body deformations for 
the normal load Fn = 3000, 5000 and 7000 N. In ah cases 
the slip s = 0.05 and the coornering angle ^ = 0. The open 
squares denote the position of rubber elements of the unde- 
formed tire body, and the filled squares underneath denote 
the position of the same tire body elements of the deformed 
tire. The maximum tire body displacements are 0.92, 1.39 
and 1.84 cm for the tire loads Fn = 3000, 5000 and 7000 N, 
respectively. For the rectangular contact area with contact 
pressure 0.3 MPa. For the rubber background temperature 
To = 80 °C and the car velocity 27 m/s. 



0.3 MPa, and with the slip 5 = and the coornering 
angle ^ = 5°. The maximum tire body displacements are 
1.20, 1.81 and 2.31 cm for the tire loads Fn = 3000, 5000 
and 7000 N, respectively. 

In Fig. [2T]we show the maximum friction coefficient, 
Mmax, of the /Xa;-slip curvc for the elliptic, rectangular 
and trapezoid footprints, as a function of the tire load. 
We show results for the contact pressures p = 0.1 MPa 
(upper three curves) and 0.3 MPa (lower three curves), 
but the tire body properties are assumed to be the same. 
Note that the efi"ective friction depends on the (average) 
tire-road footprint pressure p. When p increases, assum- 
ing unchanged size of the footprint, the friction decreases. 
This is one reason for why racer tires exhibit much larger 
friction than passenger car tires (the contact pressure in 
Fl tires is of order 0.1 MPa which is about 3 times lower 
than in passenger car tires). The increase in the friction 
as the tire-road contact pressure p decreases is mainly 
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FIG. 20: Snapshot pictures of the tire body deformations for 
the normal load Fn = 3000, 5000 and 7000 N. In ah cases the 
slip s = and the coornering angle = 5°. The open squares 
denote the position of rubber elements of the undeformed tire 
body, and the filled squares underneath denote the position 
of the same tire body elements of the deformed tire. The 
maximum tire body displacements are 1.20, 1.81 and 2.31 cm 
for the tire loads i^N = 3000, 5000 and 7000 N, respectively. 
For the elliptic contact area with contact pressure 0.3 MPa. 
For the rubber background temperature To = 80 °C and the 
car velocity 27 m/s. 
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FIG. 21: The maximum friction coefficient, /imax, of the fix- 
slip curve for elliptic, rectangular and trapezoid footprints. 
For the contact pressures p = 0.1 MPa (upper three curves) 
and 0.3 MPa (lower three curves). For the rubber background 
temperature To = 80 °C and the car velocity 27 m/s. 
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FIG. 22: The longitudinal tire stiffness Cl associated with the 
fix-slip curve for elliptic, rectangular and trapezoid footprints. 
For the contact pressures p = 0.1 MPa (upper three curves) 
and 0.3 MPa (lower three curves). For the rubber background 
temperature To = 80 °C and the car velocity 27 m/s. 
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FIG. 23: The maximum longitudinal tire body displacement 
for elliptic, rectangular and trapezoid footprints. For the con- 
tact pressures p = 0.1 MPa (lower three curves) and 0.3 MPa 
(upper three curves). For the rubber background temperature 
To = 80 °C and the car velocity 27 m/s. 



due to a decrease in the pressure in the macro asperity 
contact regions as p decreases: When the local pressure 
decreases the produced heating of the rubber decreases 
leading to a larger friction. 

In Fig. [22] we show for the same systems the longitudi- 
nal tire stiffness Cx associated with the //a;-slip curves. 
Finally, Fig. [23] gives the maximum longitudinal tire 
body displacement as a function of the tire load, for the 
same systems as in Fig. [2T} 

4,3 Dependence of the /^-slip curve on the car 
velocity 

Fig. [21] shows the the /x-slip curves for the rectangular 
footprint (20 cm x 10.2 cm) shown in Fig. [5J and for 
the car velocities Vc = 10, 30 and 40 m/s. We show 
results both for the tire load Tn — 2000 N and p = 




FIG. 24: The /i-slip curves for a rectangular footprint (20 cmx 
10.2 cm) and for the car velocities Vc = 10, 30 and 40 m/s. 
For Fn = 2000 N and p = 0.1 MPa (top three curves) and 
Fn = 6000 N and p = 0.3 MPa (lower three curves). For the 
rubber background temperature To = 80 °C. 
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FIG. 25: The /^-slip curves for the car velocity Hc = 10 and 
27 m/s. The maximum /i* of the /i-slip curve, and the slip s = 
s* where the maximum occur, depends on the car velocity Wc. 
In the present case s* « 0.057 for both velocities. The ABS 
control algorithm should increase the braking torque when 
s < s* and reduce the braking torque when s > s* . 



0.1 MPa (top three curves) and for Fn = 6000 N and 
p — 0.3 MPa (lower three curves). Note that as the 
car velocity decreases the longitudinal tire stiffness Cx 
decreases and the maximum of the /i-slip curve increases. 
This is in accordance with experimental observations for 
passenger car tires (see, e.g.. Fig. 8.65 in Ref. 3)- When 
the contact pressure and the load both increases, in such 
a way that the contact area stays constant, both the tire 
stiffness and the maximum of the /i-slip curves decreases. 

5. ABS braking simulations 

The theory developed above may be extremely useful 
to design or optimize control algorithms for ABS brak- 
ing. Here we will present results using the two simplest 



possible control algorithms. In both cases the braking 
torque is changed (increase or decrease) in steps of AM 
at time i„ = nAt {n — 1,2, ...). The first algorithm (a) 
assumes that the slip s* where the friction is maximal is 
known (and constant in time). In this case the braking 
torque is increased if the slip s(tn) at time i„ is below 
s* and otherwise it is decreased, see Fig. [25l One prob- 
lem here is that the slip s* depends on the car velocity 
which change during the braking process. However, in 
the present case s* « 0.057 nearly independent of the 
car velocity for 10 m/s < Uc < 27 m/s. 

In the second control algorithm (b) we assume that 
s* is unknown. Nevertheless, by registering if the (lon- 
gitudinal) friction Fx{t) is increasing or decreasing with 
time we can find out if we are to the left or right of the 
maximum at s = s*. That is, if 



Fxitn) > Fxitn-l) and s{tn) < s{tn-l) 



or if 



Fxitn) < Fxitn-l) and s{tn) > s{tn-l) 

then we must have s(t„) > s* and the braking torque at 
time tn is reduced, otherwise it is increased. Here Fx{tn) 
is the longitudinal friction force and s(i„) the slip at time 
t„ = nAt (n = 1,2,...). 

We now present numerical results to illustrate the two 
ABS braking algorithm. Let M be the mass-load on a 
wheel and / the moment of inertia of the wheel without 
the tire. We assume for simplicity that the suspension 
is rigid and neglect mass-load transfer. The equations 
of motion for the center of mass coordinate x{t) of the 
wheel, and for the angular rotation coordinate (j){t) are: 



Mi = Fr\ 



I(j) = Mriin - Mb 



(2) 



(3) 



where Frim is the force acting on the rim, Mb is the 
braking torque and Mj-im the torque acting on the rim 
from the tire (for constant rolling velocity F^-^^ — F[ is 
the tire-road friction force and Mj-im — RFf, where R is 
the rolling radius, but during angular accelerations these 
relation no longer holds because of tire inertia effects) We 
have used M = 360 kg and / = 0.4 kgm^. 

We assume first the control algorithm (a). We take 
AM = 200 Nm and At = 0.03 s, and we assume (see Fig. 
[25| s* « 0.05. In Fig. [26] we show (a) the car velocity 
Vc and the rolling velocity wr,, (b) the longitudinal slip 
and (c) the braking torque as a function of time. The 
time it takes (i = 1.775 s) to reduce the car velocity 
from vq = 27 to vi ~ 10 m/s correspond to a friction 
coefficient fi — {vq — vi)/gt — 0.976 which is ^ 13% 
smaller than the friction at the maximum of the /i-slip 
curve, which varies between /^max = 1.14 and 1.11 as the 
car velocity changes from 27 to 10 m/s. The slope of the 
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FIG. 26: (a) The car velocity Vc and the roUing velocity vb. as 
a function of time i. (b) the slip and (c) the braking moment 
as a function of time t. For ABS braking using algorithms a 
(see text for details). 



car-velocity line in Fig. [SGla*) for t > 0.2 s correspond 
to the friction coefficient 1.02, which is larger than the 
(average) friction calculated from the stopping time. The 
slightly smaller friction obtained from the stopping time 
reflect the (short) initial time interval necessary to build 
up the braking torque. 

In Fig. [27] we show results for the ABS control al- 
gorithm (b). Note that it takes t — 1.949 s to reduce 
the car velocity from vq — 27 m/s to vi = 10 m/s. 
This correspond to the effective friction coefficient /i = 
(vo — vi)/gt — 0.889. The maximum in the /i-slip 
curve (see Fig. [25|) depends on the car velocity but is 
about /Ltmax = 1.14 for Vc = 27 m/s and about 1.11 for 
Vc = 10 m/s so the ABS braking control algorithm used 
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FIG. 27: (a) The car velocity Vc and the rolling velocity vn as 
a function of time t. (b) the slip and (c) the braking moment 
as a function of time t. For ABS braking using algorithms b 
(see text for details). 



above could still be improved. Note also that the wheel 
tend to lock about 3 or 4 times per second. This is in 
good agreement with ABS braking systems presently in 
use. However, since the speed of cars are usually not 
known during ABS braking, braking control algorithms 
used in most cars today determine the braking torque 
only from the wheel rotation acceleration. This is possi- 
ble because, as shown in Fig. BTT a). as the wheel tend to 
lock, the rotational velocity very rapidly decreases, and it 
this point the ABS system decreases the braking torque. 
Note that the (average) of the slip in Fig. [26l b) and 
I27r b) increases with increasing time or, equivalently, de- 
creasing car velocity. This is mainly due to the fact 



that the time it takes for the wheel to lock, when the 
slip s > s*, decreases as Vc decreases. Thus, during 
the time period Ai between two changes in the brake 
torque the maximal slip (corresponding to the minimal 
rolling velocity) will increase as Vc decrease. This is easy 
to show mathematically. Since the car velocity changes 
slowly compared to the rolling velocity, from the defini- 
tion s = {vc — ur)/i'c we get 



dt 



ds 
'di 



If we approximate the /i-slip curve for s > s* with a strait 
line, 



Moff ~ Ato - A^s, 



we get from (3) 

dt^ ~ R dt 
or 



Iv^ ds 



Mg[fio - A^s] - Mb 



ds 
'di 



= -A + Bs 



where A = {MgR^o ~ Mb){R/Ivc) and B = 
A/i(Mgi?^//wc). Since A and B can be considered as 
constant during the time interval between the changes in 
the braking torque, we get 
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FIG. 28: Dynamical /^-slip curves for ABS braking using al- 
gorithms a (top) and b (bottom). The green curve is the 
steady-state /i-shp curve for the car velocity Vc — 27 m/s. 



where 



A _ 1 
B ~ A^ 

It is easy to show that 

A 



Mo 



Mb \ 
MgRj 



m -^ = im - s*] 



M* - M, 



B 



A^MgR 



where Mg = Mg{^Q~ Afis*) is the braking torque neces- 
sary in order to stay at the maximum in the /Lt-slip curve. 
If Mb < MjJ and s(0) > s* we have s(0) - A/B > and 
during the time interval Ai the slip will increase with 
[s(0) - A/B]exp{BAt). Since BAt -- l/wc the maximum 
slip will increase exponentially (until the wheel block, 
corresponding to s = 1) with the inverse of the car ve- 
locity. This behavior (i.e., the increase in the slip with 
decreasing car velocity) can be seen in Fig. [26lb) and is 
even stronger for the second ABS control algorithm [fig. 

HZIb)]. 

In Fig. [28] we show the /j,-slip curve during stationary 
slip (green curve) and the instantaneous effective fric- 
tion coefficient fJ.cs{t) = Fx{t)/F^ during braking (blue 
curve). 
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FIG. 29: The car velocity Vc as a function of time t during 
ABS braking using two algorithms a and b. The procedures a 
and b results in the time periods 1.76 and 1.95 s for reducing 
the car velocity from 27 to 10 m/s. The effective friction 
values 0.976 and 0.889 are both smaller than the maximum 
kinetic friction which is 1.098 and occur at the slip velocity 
0.0316 m/s. 



The blue and red curves in Fig. [29] show the car veloc- 
ity using the ABS control algorithms (a) and (b). It is 
clear that the ABS control algorithm (a) is more effective 
than algorithm (b), but algorithm (a) assumes that s* is 
known and remains constant during the braking event. 

The ABS braking control algorithms used today usu- 
ally assumes that only the wheel rolling velocity VB,{t) 
is known. Basically, whenever a wheel tend to lock-up, 
which manifest itself in a large (negative) wheel angu- 
lar acceleration, the braking torque is reduced. These 
ABS braking control algorithms (e.g., the Bosch algo- 
rithm) are rather complex and secret. The calculations 
presented above can be easily extended to such realis- 
tic ABS braking control algorithms and to more com- 
plex cases such as braking during load fluctuations (e.g., 
braking on uneven road surfaces) and switching between 
different road surfaces (by using different road surface 
power spectra's during an ABS braking simulation). 

6. Summary and conclusion 

In this paper we have proposed a simple rubber fric- 
tion law, which can be used, e.g., in models of tire (and 
vehicle) dynamics. The friction law gives nearly the same 
result as the full rubber friction theory of Ref. [6] , but 
is much more convenient to use in numerical studies of, 
e.g., tire dynamics, as the friction force can be calculated 
much faster. 

We have proposed a two-dimensional (2D) tire model 
which combines the rubber friction law with a simple 
mass-spring description of the tire body. The tire model 
is very flexible and can be used to calculate accurate fi- 
slip (and the self-aligning torque) curves for braking and 
cornering or combined motion (e.g., braking during cor- 
nering). We have presented numerical results which il- 
lustrate the theory. Simulations of Anti-Blocking System 
(ABS) braking was performed using two simple control 
algorithms. 
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